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Abstract. A periodic one-dimensional Schrodingcr operator is called semifinitc- 
gap if every second gap in its spectrum is eventually closed. We construct ex- 
plicit examples of semifinite-gap Schrodingcr operators in trigonometric func- 
tions by applying Darboux transformations to the Whittakcr-Hill equation. 
We give a criterion of the regularity of the corresponding potentials and inves- 
tigate the spectral properties of the new operators. 



1. Introduction 

It has been well-known since the work of Floquet and Bloch that the spectrum 
of the Schrodingcr operator 

L = -d^ +u[x) 

with a smooth real periodic potential u(x + T) = u(x) has a band structure. In 
the 1970s starting with Novikov's famous work [15] a beautiful finite-gap theory, 
linking the spectral theory with classical algebraic geometry, was developed. It 
turned out that all operators having only a finite number of gaps in their spectrum 
can be described explicitly in terms of hyperelliptic Riemann's theta functions [llj . 
The simplest one-gap operator is a particular case of the famous Lame operator 

L=-^+2 P (x), 

where p{x) is the Weierstrass elliptic function (shifted by the imaginary half-period 
to make it non-singular). None of the smooth real periodic finite-gap potentials 
(except constants) can be expressed in terms of elementary trigonometric functions. 

On the other hand, there are smooth trigonometric potentials, for which every 
second gap will be eventually closed. We call such potentials semifinite-gap. A 
trivial example can be given by any potential, which has a period T/2, but is 
considered as T-pcriodic. One can easily check that for such a potential all the odd 
gaps are closed. In fact, all the potentials with odd gaps closed are of this form 
(see Theorem XIII.91 in [17]). 

A more interesting example is given by the so-called Whittaker-Hill potential 

u(x) = A cos 2x + B cos Ax, 

with 

A = -Aas, B = -2a 2 , 

for any real a and natural s. Namely, it is known (Magnus- Winkler [14], Djakov- 
Mityagin [6]) that for odd s = 2rn + 1 all the even gaps except the first m are closed 
and for even s = 2m the same is true for odd gaps. 
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In this paper we construct new explicit trigonometric examples of semifinite-gap 
potentials by applying the Darboux transformation to the Whittaker-Hill operator. 
The main issue here is the regularity of the corresponding potentials, which we 
discuss in detail. 

We note that all non-constant trigonometric potentials, which are the results of 
Darboux transformations applied to zero potential, arc known to be singular (in 
contrast to the hyperbolic case, when we have many non-singular multisoliton po- 
tentials), so the existence of a large class of regular periodic potentials in elementary 
functions in the semifinite-gap case seems to be interesting. 

We should mention here that for the hyperbolic version of Whittaker-Hill op- 
erator with u(x) = 4ascosh2x + 2a 2 cosh Ax the Darboux transformations were 
considered by Razavy in |16j , but the general regularity question was not addressed 
in that paper (which is different in the hyperbolic case anyway). In the framework 
of quasi-exactly solvable systems the Whittaker-Hill equation has appeared as type 
X in Turbiner's list [19] and was discussed in details in [7]. 



2. Whittaker-Hill equation 
The Whittaker-Hill equation 

(1) - ip" - (4ascos2a; + 2a 2 cos Ax)ip = Xtp, 

is an eigenvalue problem for the following 7r-periodic Schrodinger operator 

(2) L = -D 2 - (Aas cos 2x + 2a 2 cos Ax) , 

where D = It depends on two real parameters a and s. Change of the sign of 
a is equivalent to the shift x — > x + ir/2 and to the change of sign of s, so we can 
assume without loss of generality that s > 0, a > 0. When a = we have the zero 
potential. 

This equation appeared in the work of Liapunov [13] and was studied extensively 
by Whittaker [21] (see also Ince [10]). Here we follow the terminology of the book 
[14] by Magnus and Winkler. The material of this section is mainly based on this 
book and on the recent paper [6] by Djakov and Mityagin. 

The substitution 

(3) ^ = ipe acos2x 
carries (fTJ) into 

(4) ip" - Aasm2xtp' + [\ + 2a 2 +A{s -l)a cos 2x]ip = 0. 
It can be rewritten as 

(5) Kip = v(p, 
where 

(6) K = -D 2 + 4a(sin 2x)D - A(s 

and v = A + 2a 2 . 

Consider 7r-periodic solutions to ([4]) of the form 

oo 

(7) A) + ]T a 2 cos2na;, ip odd 

n=l 



— l)acos2x 



oo 

£?2n sin 2nx. 

n=l 
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Substituting Q into © we end up with the following recurrence relations 

(8) vA Q + 2a(s + l)A 2 = 

(9) 4a(s - l)A + {y - A)A 2 + 2a{s + 3)A 4 = 

(10) 2a{s -2k + l)A 2k - 2 + Ak 2 )A 2k + 2a{s + 2k + l)A 2k+2 = {k > 2) 



in the even case and 
(11) 



{u - 4)B 2 + 2a(s + 3)S 4 = 



(12) 2a(s -2k + l)B 2k -2 + (v - Ak 2 )B 2k + 2a(s + 2k + l)B 2k+2 = (k > 2) 
in the odd case. Similarly, for anti-periodic solutions 



(13) 



^ C 2n +i cos(2n + l)x, Lp dd = D 2n+ i sin(2n + l)a 



we have the following recurrence relations: 

(14) (v + 2as- l)Ci + 2a(s + 2)C 3 = 

(15) 2a(s - 2k + 2)C 2fe _ 3 + {v- (2k - l) 2 )^-! + 2a{s + 2k)C 2k+1 =0 (k > 2) 
and 

(16) (v - 2as - l)Di + 2a{s + 2)D 3 = 

(17) 2a{s-2k + 2)D 2k - 3 + {v-{2k-lj 2 )D 2k - l + 2a{s + 2k)D 2k+1 =0 (k > 2). 

We will start first with the periodic case. The following proposition [Til [5] 
explains the special role of the integer parameters s in the theory of the Whittaker- 
Hill equation. 

Proposition 1. Suppose that for some A the equation ^ has two non-zero solu- 
tions Ffy with A 2n ,B 2n decaying exponentially fast as n — ► oo. Then s must be an 
odd integer. 

Proof. We have from (JTUJ) and (0 for k > 2 

a k A 2k - 2 + b k A 2k + c k A 2k+2 = 0, 
a k B 2k ^ 2 + b k B 2k + c k B 2k+2 = 0, 

where a k = 2a(s — 2k + 1), b k = v — k 2 and c k = 2a(s + 2k + 1). By multiplying 
the first equation by B 2k and the second by A 2k and subtracting, we have 



a k 



A 2k - 2 A 2k 
B 2k - 2 B 2k 



Ck 



A 2k A 2k+2 
B 2k B 2k+2 



0. 



Introducing A k 



A 2k - 2 A 2k 
B 2k - 2 B 2k 



we have 
a k A k = CfcA fc+ i 
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and hence 



A,, 



k+r+l 



■ Afe+ r +i 

flfcOfc+l • ■ ■ &k+r 

(s + 2k + l)(s + 2k + 3) . . . (s + 2k + 2r + 1) 
(s -2k + l)(s - 2fc - 1) . . . (s - 2k - 2r + 1) 

2 r+1 {k + s/2 + l/2)(fc + s/2 + 3/2) . . . (k + s/2 + r + 1/2) 
(-2) r+1 {k - s/2 - 1/2) (fc - s/2 + 1/2) . . . (k - s/2 + r - 1/2) 
Nr+1 T(fc + s/2 + r + 3/2)r(fc - g/2 - 1/2) 



A 



fc+r+l 



T(fc + s/2 + l/2)r(fc - s/2 + r + 1/2) 

where T denotes the classical Gamma function [22] and we have used the functional 
relation 

T(z + r + 1) 



T(z) 



z(z + l)...(z + r). 



We rewrite this as 



r(fc + s/2+l/2) _ r+i r(fc + s/2 + r + 3/2) 
r(fc-s/2-l/2) fe_l J r(fc-s/2 + r + l/2) fc+r+1 ' 



Now from Stirling's formula [22] we have the asymptotic relation 

1 



lim 

t — >+oc 



r(t- P ) 



for any fixed p. Since Ai n and Bi n decay exponentially fast, so does A„. By setting 
t = k + r + 1, p = s/2 + 1/2, and letting r — > oo in ([2]) we see that the right hand 
side goes to zero, which implies that the left hand side (which is independent of r) 
must be zero: 

r(fc + s/2+l/2) 
r(fc-s/2-l/2) k 

or, more explicitly 

(fc - s/2 - 1/2) (fc - s/2 + 1/2) . . . (fc + s/2 - 1/2) (fc + s/2 + l/2)A fe = 

Hence, we either have A& = for all fc > 2, or one of the brackets must vanish. In 
the first case it is easy to see that either ip even or (p dd must be identically equal to 
zero. In the second case s is an odd integer. □ 



Let now s = 2m + 1 be a positive odd integer. In that case a m +i = 2a(s — 
2rn — 1) = 0, so the infinite tri-diagonal matrices corresponding to the relations 
(|10p and (|12| become reducible. Introduce the corresponding finite-dimensional 
parts, corresponding to the terminating sequences A^n = i?2n = for n > m : 

\ 

... 
c 2 ... 

n — 1 ^m— 1 C-m — 1 

a m b m J 



(18) 



K° = 



(bo 


CO 





Ol 


h 


Cl 





a 2 


62 
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and 





/hi 


Cl 


ii 


\ 








&2 


C2 


... 









«:S 


6 3 


c 3 ... 




(19) Kl = 






















^m — 1 b rn —\ C m _i 












a m 6 m y 




where bk = v — 4fc 2 , Ck = 


2a(s 


+ 2k 


+ 1) 


for all fc, = 2a(s - 


- 2fc + 


and ax = 4a(s — 1). 










Let 












,a) = 


detiC 


^(".a) = deti^ 


and 










So (a) = {v : 




a) 


0}. 


Si (a) = : <^(f,a 


) = 0} 



be the corresponding eigenvalues. 

Since a > all the off-diagonal entries of the matrices and are positive. 
From the general theory of tri-diagonal (Jacobi) matrices (see e.g. [TB]) it follows 
that all the eigenvalues v are real and distinct. Moreover, these sets So and Si are 
interlacing: between any two zeroes of So there is exactly one zero of Si . Although 
this fact is well-known for the convenience of the reader we give the proof in the 
Appendix 1. 

Remark. In modern terminology S = So U Si corresponds to the solvable part of 
the spectrum of the Whittakcr-Hill operator, which in the case when s is an integer 
belongs to the class of quasi- exactly solvable operators [19(120] . In the theory of such 
operators the polynomials are sometimes called Bender-Dunne polynomials [2] • 

A crucial observation due to Magnus- Winkler [14] (see also [6]) is that the solv- 
able part coincides with the simple part of the periodic spectrum. More precisely, 
we have the following 

Proposition 2. Let s = 2m + 1, m € Z+. If v belongs to the periodic spectrum of 
the Whittaker-Hill equation and has multiplicity 1, then c£S. 

Proof. Let us assume that v (jt S and that we have a non-zero even periodic solution 
of the form 



An 



fe=l 



Ank cos 2kx. 



Since v ^ S the coefficient A 2m+2 7^ (otherwise, we have A 2 k = for all fc > to 
and thus v € S ). We now construct a second, odd periodic solution for the 
same v. First we set B 2n = A 2n for n > m, which is OK since the recurrence 
relations (fTu} and (fl"2"|) agree for k > 2. Now, since v ^ Si the m x to matrix Kq 
is invertible. This means that we can reconstruct uniquely the beginning of the 
sequence B = (B 2 ,B4, . . . , B 2m ) as 



/ 







(20) 



B = {Kl 



(i 
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This means that we have independent periodic solutions with the same v 1 so v has 
multiplicity 2. Contradiction means that v G So in that case. A similar argument 
in the case of odd periodic solution shows that v G S\. □ 

Similar results are true for even s = 2m and anti-periodic spectrum. From 
(fl4|) .(fT ^ . (fT6|) . (|17|) the corresponding eigenvalues coincide with the eigenvalues of 
the following matrices: 



(21) 



fbf 


Cl 







... \ 




h 


C2 











0,3 


&3 


C3 













0-771— 1 


^m— 1 Cm-1 


\o 












where ct = 2a(s + 2k), at = 2a(s — 2fc + 2)for all k, bk 



(2fc-l) 2 for k > 2 and 



1 ± 2as. The corresponding sets are also interlacing (see Appendix 



1). We denote all the eigenvalues from the set S = S + U S in increasing order as 

V\ , . . . , V<x m ■ 

The degeneracy of the eigenvalues v ^ S is called the coexistence property in 
|14j . This leads to the following spectral property of the Whittaker-Hill operator. 
For the general theory of periodic Schrodinger operators we refer to the classical 
Reed-Simon book [T7] . 

Let us choose a basis ipi( x )i ^2 (a?) of the solution space of the Whittaker-Hill 
equation (JTJ) and define the corresponding monodromy matrix by 

' tftlW(A)( tfl 

Taking wronskians of both sides we see that detM(A) = 1. The trace of the mon- 
odromy matrix A(A) = trM(X) is independent of the choice of the solutions ipi, ip 2 
and sometimes is called Hill's discriminant (in Russian literature the term Liapunov 
function is also used). The Floquet multipliers /ii,/i2 are the eigenvalues of M(A), 
which satisfy the characteristic equation 

/1 2 - A(A)/i + 1 = 0. 

The continuous spectrum of L consists of the segments of A G K such that corre- 
sponding solutions are bounded on the real line. This happens when \fii\ = 1, or 
equivalcntly 

|A(A)| <2. 

The complement is the union of intervals called gaps. They correspond to the part 
of graph of A(A) lying outside of the spectral strip — 2 < A < 2. 

In general (for example, for the Mathieu operator with u = A cos 2x) this graph 
intersects each line A = ±2 during each wave of oscillation (see Fig. 1, left). In 
the finite-gap case starting from some point the graph starts to touch these lines 
but not intersect them (see the one-gap case u(x) — 2p(x) on Fig.l, right). 

In the scmifinitc-gap case this happens only for one of these lines. For example, 
for the Whittaker-Hill operator the graph eventually touches the periodic line A = 2 
for odd s and the anti-periodic line A = —2 for even s (see Fig. 2). The case of small 
a (see below) shows that the open gaps are precisely the first ones; by continuity 
arguments this is true for all real a. 
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Figure 1. Graphs of the Hill's discriminant. Left: generic peri- 
odic case. Right: one gap case 




Figure 2. Hill's discriminant in the scmifmitc-gap case. Left: 
s = 2. Right: s = 3. 

We summarize the results of this section as follows. 

Theorem 1. 14, 6] The Whittaker-Hill operator has all even gaps closed except 
the first m when s = 2m + 1 and all odd gaps closed except the first m when s = 2m. 

Now we are going to show that there are many other semifinite-gap operators 
with potentials which can be expressed in terms of trigonometric functions. The 
idea is standard: to apply the Darboux transformation [5]. The main problem is 
how to get non-singular potentials. We are going to see first what happens for small 
s and a. 

3. Examples and the small a case 

In the first non-trivial case s = 1 the corresponding operator ^ simplifies to 

K = -D 2 +4asm2xD 

and clearly has an eigenfunction ipa = 1 with the eigenvalue v = 0. The correspond- 
ing Whittaker-Hill operator 

L = -D 2 - (4a cos 2x + 2a 2 cos Ax) 

has an explicit ground state 

/ a cos 2x 

ipo = e 

with the eigenvalue A = —2a 2 . 
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When s = 2 we have 

K = -D 2 + 4a sin 2xD - 4a cos 2x. 

One can easily check that it has ipi = cos a; and (f2 = sins as the eigenfunctions 
with v\ = 1 — 4a and vi = 1 + 4a respectively. 
In the s = 3 case we have the matrix 



K° = 



v 8a 
8a v - 4 



so we have 



5° = v{v-4) -64a 2 , <5 x = ^-4 



and the following table 





Eigenvalue 




Eigenfunction 


^0 


2(1 - VI + 16a 2 ) 


Po 


1 + ^+\^ z -i cos 2a; 


fl 


4 


Pi 


sin 2a; 


^2 


2(1 + Vl + 16a 2 ) 




l-Vi+16^ +cos2a; 



For small a we have 



vq « —16a , </?o ~ 1 + 2a cos 2a;, 
i>l =4, y?i = sin 2x, 



V2 ~ 4 + 16a , (^0 



-2a + cos 2x. 



When s — 4 the matrices arc 



j/-l±8a 12a 
4a ^ — 9 



The eigenvalues are the zeros v = v^ 2 of <5 ± = dct : 

8 ± {v) = (v-l± 8a)(u - 9) - 48a 2 = 
The corresponding eigenfunctions are respectively 

Pi 2 = i v i2 ~ ^) cos x — cos 3a;, 



0. 



When a 



Pi 2 = (^i 2 — 9) sin a; — 4a sin 3a;. 
we have the first s periodic levels of the corresponding operator 





= 0, 




= 1 


"1 




Pi 


= sin 2.t 




= 4, 


Y2 


= cos 2a; 



i/ 2ra -i = 4m , tpim-i = sin2ma3 
V"Xm = 4m 2 , if 2m = cos 2ma; 
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when s = 2m + 1 and 

i/i = l, tpi = sin a; 
Vi = 1, <^2 = cosx 

"2 m -i = (2m - l) 2 , ^2m-i = sin(2m - l)x 
U2 m = (2m — 1) 2 , <p2m = cos(2m — l)x 

for s = 2m. 

Note that since we denote the eigenvalues in increasing order it follows from 
the interlacing property that for odd s the even eigenvalues correspond to the 
even eigenfunctions and odd eigenvalues to odd cigenfunctions. For even s the 
opposite is true: odd eigenvalues correspond to even eigenfunctions and vice versa 
(see Appendix 1). 

We are going to start with this limiting case to see the properties of the wron- 
skians of the corresponding cigenfunctions. 

4. DARBOUX TRANSFORMED WHITTAKER-HlLL OPERATORS 

We recall the following classical construction going back to Darboux [5]. Let 
tf>i , . . . , ip n be a set of eigenfunctions of the Schrodinger operator 

L = -D 2 + u(x) 

with some eigenvalues k\, . . . , K n : 



-l/)'! + u(x)l/)j = Kjlpj, j = 1, 



. n. 



v 3 T uy^jyj — ruji/zj, 

The corresponding Darboux transformation of L is the operator L = —D 2 + u(x), 
where 

(22) u = u-2D 2 logW(4>i,...^ n ), 

where W(if>i, . . . , if) n ) is the wronskian of the functions ipx, . . . , ip n . The key property 
of the new operator is that the generic solutions of the corresponding Schrodinger 
equation 

can be explicitly given in terms of the solutions of the initial equation 

-ip" + u(x)tp = 

by the Crum formula [3]: 

(23) ^ = • 

This works well for all A ^ Kj, j = 1, . . .n. 

We will need also the following information about eigenfunctions for A = Kj. 
When n = 1 the inverse 

(24) & = -L 

is a solution of the transformed equation with the eigenvalue A = «i . For n = 2 one 
can check that 

(25) i'i = TTTi-, 7-r, tp2 
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satisfy the transformed equation with A = Ki and A = k>2 respectively. 
Consider now the Whittaker-Hill operator with 

u = — (Aas cos 2x + 2a 2 cos Ax) 

with odd s = 2m + 1 and the corresponding set of the explicit eigenfunctions 
^0) • • • i ip2m from the previous section. Let I = {ii ,...,£&} be the subset of 
{1, . . . , 2m} and 

W I = W^ il ,...,i> ik ) 

be the corresponding wronskian. Let us introduce also the wronskian of the corre- 
sponding functions ipj = ipje~ acos2x : 

Vj{x) = Wfa,. ..,<p ih ) = e - kacos2x W!. 

We are looking for the subsets I such that the corresponding Darboux transfor- 
mations 

(26) u = u - 2D 2 log Wi = Aa(2k - s) cos 2x - 2a 2 cos Ax - 2D 2 log Vj 

are non-singular. This clearly happens if and only if the wronskian Wj has no real 
zeros on the whole line. 

We claim that this holds for the following subsets. First of all we call a cluster 
any pair {2k — 1, 2k} as well as a single element set {0}. By definition, the cluster 
subsets I C {0, . . . , 2m} are those which consist of several clusters. In other words, 
the subset I is called a cluster when it contains an element 2k — 1 if and only if 
it contains 2k for all k = 1, . . . , m. This applies also to the case of even s = 2m, 
except that in that case there is no single element clusters since I C {1, . . . , 2m}. 

Theorem 2. For any cluster subset I the corresponding wronskian Wj(x) has no 
zeros on the real line. 

Proof. First of all this is obviously true for a = (and hence, for small a). Indeed, 
for I = {2k — 1, 2k} we have in this case 

sin 2kx cos 2kx 

1 (2k)cos2kx (-2fc)sin2fcr 

= (-2k) [sin 2 2kx + cos 2 2kx] = -2k, 

which is a constant, hence clearly non-zero for positive k. The same can be shown 
for the general cluster case. 

Now we claim that the cluster wronskian Wj remains non-vanishing anywhere 
for all (not necessarily small) real values of a. We will prove this by induction on 
the number k of elements in /. 

If k = 1 then / = {0} and Wj = ipo- But i/jq is the ground state and therefore it 
has no zeroes by a general theorem (see e.g. [17]). 

If k = 2 then I = {2k - 1, 2k} and Wi = W(ip 2 k-i,ip2k). We know that Wi(x; a) 
is a smooth function of both x and a and that it has no zeros for small a. Suppose 
that there is a positive value of a for which Wj has at least one zero and consider 
the minimal value a* of a for which is true. By definition, for a = a* there exists 
real x = Xq such that Wi(xo) = and Wj(xo) = 0. We claim that this is impossible. 
Indeed, we have 

Wj = (lp2k-1^2k ~ -02fc-i/4fc-l)' = i>2k-X^2 k - V>2feV4'fc-l = (A2fc-1 ~ ^2k)^2k-1^2k- 
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We see that Wj(xo) = implies that either tp2k-i(%o) = or ip2k(xo) = 0. Assume 
without loss of generality that ip2k(xo) = 0. Then from Wi(xo) = we must have 
either V^fcO^o) = or ip2k-i{xo) = 0. The first case is impossible since then ip2k 
must be identically zero as a solution of second order linear differential with zero 
initial data at xq. Suppose that i/^fc-i {%o) — 0, which means that ip2k-i an d ip2k 
have a common zero xq. We claim that this is impossible because of the following 
classical result. 

Sturm's Theorem [10] . Let u(x), v(x) be smooth functions on the interval 
[a, b] such that 

u"(x) = Gi(x)u(x) 
v"{x) = G 2 (x)v(x) 

with G± > G2 on [a, b]. Then between any two consecutive zeros of u there exists a 
zero of v. 

Applying this to the Whittakcr-Hill equation ip" = (u(x) — X)ip we see that since 
^2k—i < ^2k we must have a zero of ip2k-i between any two consecutive zeros of ip2k- 
We know that for a = both ip2k = cos2/ca; and ip2k-i = sin2fca; have 2k zeros on 
[0, 7r). When a is increasing from zero ip2k and ip2k-i must maintain the number 
of zeros since otherwise there must be a value of a for which the corresponding 
solution ip has a multiple zero, which is impossible for a second order differential 
equation. But this, in combination with the coincidence of the zeros of ip2k and 
ip2k-i, contradicts the interlacing property guaranteed by Sturm's theorem. This 
proves the one cluster case. 

Suppose now we have cluster / = {2k — 1, 2k, 21 — 1, 21}. Then we can do the 
corresponding Darboux transformation in two steps. First we apply it to the cluster 
Ik = {2k — 1,2k}. We know already that the result is a non-singular periodic 
Schrodinger operator having the transformed eigenfunctions ip2l-i>ip2l- Now we 
repeat the same arguments for these eigenfunctions to show that W(ip2i-i,ip2i) has 
no real zeros. But 



W($2l-\-,i>2l) 



W(lp2k-l,1p2k) 



so this implies that Wi = W(ip2k-it i>2k) ip2l-i> ^21) has no real zeros as well. Con- 
tinuing in this way we prove the theorem. □ 

Thus for s = 2m + 1 we have constructed 2 m+1 — 1 new potentials corresponding 
to the 2 m+1 — 1 cluster sets. However, some of them are equivalent. For example, 
the largest cluster set I = {0, . . . , 2m} corresponds to the potential 

L* = -D 2 + Aas cos 2x - 2a 2 cos 4x, 

which is just the same potential ([2|) shifted by x — ► x + tt/2. Indeed, in that case 

Vi = W(<po , . . . , ip2m) = W(1, sin 2x, cos 2x, . . . , sin277ix, cos 2mx) = const 

since the linear span of ipo, ■ ■ ■ , tp2m coincides with the whole space of trigonometric 
polynomials up to degree 2m. This is similar to the duality in the sextic case 
considered in [8]. 

This allows us to consider only the cluster subsets S C {1,2,..., 2m} in both 
s = 2m + 1 and s = 2m cases. The corresponding Schrodinger operators have the 
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same continuos spectrum, but different auxiliary eigenvalues tj, corresponding to 
the Dirichlct eigenvalue problem 

i'(0) = ip(ir) = 0. 

Since all our potentials are even these eigenvalues coincide with one of the ends of 
the corresponding gap. For the initial Whittakcr-Hill operator ((2]) with s = 2m + 1 
the first even 

J2i = A 2 i_l, 1 = 1, ... ,771 

coincides with the left end of the gap, since we know that the corresponding eigen- 
function is odd. When s = 2m the first odd 

72i-l = A 2l , i = 1, ... ,771 

coincide with the right end of the gaps since in that case ip2i are odd. Note that in 
the limit s — > 00, a — ► such that sa — > A/ 4 we have the Mathieu operator with 
u(x) = — A cos 2a;, so this gives us the position of all 7, in the gaps in that classical 
case. 

When we apply the Darboux transformation to the cluster {2i — 1, 2i}, as one 
can see from the formulas ([23|) and ([25| . that the parity of all eigenfunctions is 
preserved except the corresponding ip2i—x, ij>2i- This means that the corresponding 
7i switched the side of the gap, while all other remain in the same position. 

Thus we have 

Theorem 3. For any integer s and any cluster k-element subset Id {1, . . . , 2tti}, m = 
[s/2] the Schrodinger operator 

Li = -D 2 + 4a(2fc-s)cos2z-2a 2 cos4:z;--2D 2 logV7(a;) 

is a non-singular periodic operator, having the same semifinite-gap spectrum as the 
Whittaker-Hill operator (0), but different positions of the Dirichlet eigenvalues in 
the gaps determined by the set I. 

Thus for any integer s we have constructed 2 m — 1, m — [s/2] new non-singular 
isospectral deformations of Whittaker-Hill equations with potentials expressible in 
terms of trigonometric functions. Some examples and graphs of the new potentials 
are given in Appendix 2. 

A remarkable fact is that cluster sets realize all possible combinations of the 
positions of the Dirichlct spectrum in the open gaps, so this gives all corresponding 
even semifinite-gap operators. This implies that the converse to the Theorem [2] is 
also true. 

Corollary 1. IfWj(x) has no real zeros then I must be a cluster subset. 

Indeed, the spectral data (including the Dirichlet eigenvalues) for the corre- 
sponding potentials must coincide with one of the described above, and hence by 
the uniqueness theorem / should be one of the cluster sets. 

5. Appendix 1: interlacing property of eigenvalues 

Recall that n x n matrix is called a Jacobi matrix if it is tri-diagonal with pos- 
itive off-diagonal elements. Such matrices play an important role in the theory of 
orthogonal polynomials [18] . 
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Let J n be an n x n Jacobi matrix 



(27) 



Jn = 



(b 


Co 







\ 


ai 


bi 


Cl 











a 2 


b 2 


c 2 













a n -2 f>n-2 


Cn-2 







On-i b n -ij 



and J rl _i be its (n — 1) x (n — 1) submatrix 



(28) 



Jn-l 



(b 


CO 







... \ 


ai 


b x 


Cl 











a 2 


b 2 


C2 










a 


n — 


3 b n -3 C„_3 


\o 









a n -2 bn-2/ 




Of Jn 


anrf J„. 


-i 


are simple and i 



Proof. The proof of the simplicity of spectrum of Jacobi matrices is well-known: it 
follows immediately from the tri-diagonal form of the linear system ( J n — XI n )v = 
that for given eigenvalue A the corresponding eigenspace is one-dimensional. 

Then by Afc(A) we denote the characteristic polynomial for the corresponding 
k x k submatrix Jk : 

A fc (A) = det(A7 fc - J k ). 

By definition we take Ao = 1. We are going to prove that the zeroes of A^+i and 
Afc are interlacing. More precisely, we will show that if Ai < X 2 < . . . < \k+i be 
the zeros of Afc+i(A) then each interval [Ai, A,+i], i = 1,2, ... ,k contains exactly 
one zero of Afe(A). 

We have the recurrence relation 

A fe (A) = det(\I k - J k ) = (A - 6 fc _!)A fc _i(A) - c fc _ iafe _ 2 A fe _ 2 (A) for all k > 2, 

which is of the form 

(29) A fc (A) = (A k \ + B^Ak-^X) ~ C k A k - 2 (\) 

with Ak = 1 > 0, B k = —bk-i and C k = c k -\a k - 2 > since J k is Jacobi. We claim 
that the following identity holds 



(30) 

where 7; = C l+2 C i+ 3 
relation (1291): 



7oA (a;)A (y) + 7l A 1 (x)A 1 ( 2 /) + . . . + -f k A k (x)A k (y) 
= A k+1 (x)A k (y) - A k (x)A k+1 (y) 
x-y 

. Cfc+i and jk = 1- This follows directly from the recurrence 



A k+1 (x)A k (y) - A k (x)A k+1 (y) 
= {(Ak + ix + Bk+i)Ak(x)-Ck+iAk-i(x)}Ak(y) 

-A k (x){(A k+1 y + B k+l )A k {y) - C k+ iA k -i{y)} 
= A k (x - y)A k (x)A k (y) + C k+ i{Ak(x)A k -i{y) - A k -i(x)A k (y)} 
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Recalling that A k = 1, this becomes 

A k+1 (x)A k (y) - A k (x)A k+1 {y) 



x-y 
A k (x)A k (y) + C k+1 



{A fc (x)A fc _!(y) - A fc _!(x)A fc (y)} 
x-y 

By induction, we obtain (|30|) . We notice that, for the special case x = y, we have 

7o{A (x)} 2 + 7 i{Ai(z)} 2 + • ■ • + lk{A k (x)} 2 
= A' k+1 (x)A k (x) - A' k (x)A k+1 (x) 

Since all C 2 : > 0, it follows that for all j we have 7j > 0. Hence, 

(31) A' k+1 (x)A k (x) - A' k (x)A k+1 (x) > 

because Aq(x) = 1 and is therefore positive. 

Now, if £ and r\ are two consecutive zeros of A k [x), we have A' k (£)A' k (r)) < 0. 
On the other hand, (J3IJ) yields -A' fc (£)A fe+ i(£) > 0, -A' k (n)A k+1 (n) > 0, so that 
Afc + i(£)Afc + i(?7) < 0. This indicates an odd number of zeros of A k+ i(x) in the 
interval £ < x < r\. In particular, there is at least one. Now let £ = x k be 
the greatest zero of A k (x); then A' k (£_) > 0, and (f5Tj) yields A k+ i(£) < 0. Since 
Afc+i(&) is positive, we obtain at least one zero of A k+ i(x) on the right of £ = x k , 
and similarly at least one on the left of the least zero x\ of A k (x). Consequently, 
we can only have one zero of A k -\-i(x) between x v and x^+i, v — 1,2, ... ,k. By 
interchanging the roles of A k (x) and A k+ i(x), we can prove that there exists one, 
and only one, zero of A k (x) between every zero of A k+ i(x). □ 

Applying this result to the matrix (fTB]) we have as a corollary that the periodic 
eigenvalues of Whittakcr-Hill operator with s = 2m + 1 from the sets So and Si 
are interlacing. 

To prove a similar result for even s and anti-periodic spectrum we should modify 
the arguments. Note that the corresponding matrices have the same size mxm 
and differ only in the first element (see (JH])). By the Weinstein-Aronszajn formula 
[T2] their characteristic polynomials A ± = dct(A/ r „ — K^) are related by 



A+ = A" 



4asAi, 



where Ai is the characteristic polynomial of the corresponding submatrix 
(hi c 2 \ 



(32) 



a 3 6 3 c 3 ... 
... 

— 1 — 1 Cm — 1 

\0 a m b m J 



By theorem 3] the zeros of A and Ai are interlacing, and hence the same is true 
for the zeros of A - and A + . Thus we have 

Theorem 5. The 'periodic eigenvalues of the Whittaker-Hill operator with s = 
2m + 1 from the sets So and S± are interlacing. The same is true for the sets S + 
and S~ for even s and anti-periodic eigenf unctions. 
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Note that since a > the set S~ is shifted to the right compared to S + , so 
the eigenvalues corresponding to the even eigenfunctions are smaller than their odd 
counterparts. 



6. Appendix 2: Examples and graphs 

We give first the explicit form of the Darboux transformed Whittakcr-Hill po- 
tentials 

Vi = 4a(2fc - s) cos 2x - 2a 2 cos Ax - 2D 2 log Vj{x) 

in the simplest cases. 

The first new example comes from s = 3 and the cluster set I = {0}. It corre- 
sponds to the ground state 



i'o = (1 + C(a)cos2x)e acos ^, C{a) 



s/l + 16a 2 - 1 



4a 



with eigenvalue v$ = 2(1 — \/l + 16a 2 ). The transformed potential has the following 
form 

v = -4a cos 2x - 2a 2 cos4x - 2D 2 log(l + C(a) cos2x), 
or, more explicitly, 

8C(a)(C(a) +cos2x) 



(33) vq = — Aa cos 2x — 2a cos Ax + 



(1 + C(a) cos2x) 2 



The corresponding Schrodinger operator is periodic, non-singular and has the ground 
state 

g — ct cos 2x 

^° = 1 + C(a) cos2x' 

When s — 4 we can take cluster / = {1,2}. The corresponding anti-periodic 
eigenfunctions are 



s • n a Vl - 2a + 4a 2 + a - 1 

tpi = sin x + A{ a) sm 6x, A= , 

3a 



y/l + 2a + 4a 2 - a - 1 

ifi2 = cos x + B(a) cos 6x, a = . 

3a 

The corresponding wronskian (up to a sign) is 

Vi = l + 3AB + 2(A + B)cos2x + (B-A) cosAx 
and the new potential is 

(34) v 12 = -2a 2 cos Ax - 2D 2 log[l + SAB + 2(A + B) cos 2x + (B - A) cos Ax] . 

The corresponding graphs as well as some of the transformed potentials in the 
case when s = 5 are shown in Fig. 3 and 4. These were created using a programme 
written in Maple. 
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Figure 3. New scmifinitc-gap potentials. Left: s 
{0}, a = 1. Right: s = 4, i" = {1, 2}, a = 1. 




Figure 4. New semifmite-gap potentials. Left: s 
{3, 4}, a = 1. Right: s = 5, J = {0, 3, 4}, a = 1. 



= 5,7 



7. Concluding remarks. 

This paper was motivated by our attempts to study the integrability of the 
quantum Neumann system. In the simplest case we just have the quantum pen- 
dulum described by the Mathieu equation with u(x) = Acos2x. It is well-known 
that this equation can be considered as a limit of the Lame equation with potential 
u(x) = s(s + l)p(x), when the parameter s goes to infinity and the elliptic curve de- 
generates to the rational one (see Section 15.5 in [3]). It seems to be surprising that 
the elliptic deformation turned out to be simpler than the original trigonometric 
version. 

What we have shown in this paper (and what had actually been discovered a long 
time ago by Whittaker [21] ) is that one need not consider elliptic generalisations 
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but instead simply add the second harmonic to the potential in order to see some 
integrability. Indeed, one can consider the Mathieu case as a limit of the Whittaker- 
Hill equation ((T|) when s — > oo, a — * such that as — > A. 

An interesting feature of the Whittaker-Hill case is the appearance of the hier- 
archy of the orthogonal trigonometric polynomials. Indeed, from the orthogonality 
of the eigenfunctions of the Whittaker-Hill operator ([2]) it follows for every odd 
s = 2m + 1 that the corresponding eigenfunctions ipk{x), k = 0, . . . , 2m of the op- 
erator ^ form a basis in the s-dimensional space V m of trigonometric polynomials 
of degree < m, which is orthogonal with respect to the scalar product 

(35) (<Pi,<P2)= l\i{x)^ 2 (x)e 2acos2x dx. 

Jo 

After the change of coordinate z = cos 2x the even eigenfunctions tp2i (x) = pi (cos 2x) 
become the usual polynomials, which are orthogonal with respect to the measure 

(36) dft{z) = -j== . 

However, in contrast to the usual case |18j they have the same degree m if a =/= 0. 
When a = we have the classical Chebyshev polynomials: tp 2 i{x) = cos2Zx = 
Ti (cos 2x). The usual orthogonal polynomials with respect to the measure (|36|) 
have been recently studied by Basor, Chen and Ehrhardt in [1]. Some non-standard 
orthogonal polynomials arc discussed by Gomez-Ullate, Kamran and Milson in the 
recent paper [5]. 

The situation here is similar to the paper [8], where another quasi-exactly solv- 
able case of sextic potential u = x 6 — (m + l)x 2 was considered. The corresponding 
Darboux transformed operators have trivial monodromy in the complex domain and 
interesting geometry of the singular sets. It is worthy to look at these properties 
in the Whittaker-Hill case as well. 
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